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Abstract 

We consider a class of stochastic processes modeling binary interactions in an 
A^-particle system. Examples of such systems can be found in the modeling of 
biological swarms. They lead to the definition of a class of master equations that 
we call pair interaction driven master equations. We prove a propagation of chaos 
result for this class of master equations which generalizes Mark Kac's well know 
result for the Kac model in kinetic theory. We use this result to study kinetic limits 
for two biological swarm models. We show that propagation of chaos may be lost 
at large times and we exhibit an example where the invariant density is not chaotic. 
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1 Introduction 

This paper is devoted to the passage from stochastic particle systems to kinetic equations 
when the number of particles tends to infinity. We are specifically interested in pair- 
interaction processes that are inspired from biological swarm models. We start from the 
level of the master equation which describes the evolution of the iV-particle probability 
distribution of the system. The master equation is posed on a large dimensional space 
consisting of an iV-fold copy of the phase space. By contrast, the kinetic equation provides 
a reduced description based on the single particle distribution function posed on a single 
copy of the phase space. To show that this reduced description is valid, one needs to show 
that the particles become statistically independent in the limit N ^ oo; that is, that the 
iV-particle probability distribution becomes a iV-fold tensor product of the single particle 
distribution function. More precisely, the result to be shown is that, if the particles are 
initially pairwise independent, the time evolution approximately propagates this pairwise 
independence, and does so exactly in the large limit; this is called "chaos propagation". 
This property of an A^-particle stochastic evolution is crucial to the existence of a kinetic 
description. 

Kinetic models derived from particle systems abound in the literature. However, only 
in very few cases has the propagation of chaos been proved, and hence only in very 
few cases have these models been mathematically derived from an underlying particle 
dynamics. The most emblematic kinetic model, the Boltzmann equation has received 
most of the attention. Following seminal works by Kac |26l |27] and McKean |33], the 
first rigorous establishment of the Boltzmann equation is due to Lanford [2^ ISUj \3T\ 
and King [2S] for Hard- Sphere dynamics and hard potentials (see also [231 1211 SO] and 
[T9| [20| |2T| |3ll E]). A new approach yielding global- in-time results has been recently 
developed by Mischler, Mouhot and Wennberg [361 ETJ [38] • Kac proposed a caricature of 
the Boltzmann equation leading to the Kac kinetic equation [26]. Propagation of chaos 
for the Kac master equation and the related question of gap estimates have received a 
great deal of attention (see e.g. [3 El El EDI [IHl 123 132]) • 

In this paper, our goal is to investigate a class of processes which are inspired from 
biological swarm models. A first example is the BDG model, named after Bertin, Droz 
and Gregoire [2]. This model is intended to be the kinetic counterpart of the Vicsek 
particle system [13]. In the Vicsek model, particles moving with constant speed update 
their velocity by trying to align with the average velocity of their neighbors. In [2], the 
authors propose a kinetic formulation of a binary collision process which mimics this 
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alignment tendency: at each collision, the two particles change their velocity to their 
average velocity up to a certain noise. One of the goals of the present paper is to provide 
a rigorous justification of this kinetic model, at least in the space-homogeneous case. 

Here, we also propose a different and to our knowledge original binary collision process 
which mimics the Vicsek alignment dynamics. In this process called 'Choose the Leader 
(CL)', one of the two colliding particles (the follower) decides to take the velocity of the 
other one (the leader) up to some noise. The choice of the leader and the follower is 
random with equal probabilities. We propose a kinetic formulation of the CL process and 
rigorously establish it in the space homogeneous case. One of the advantages of the CL 
process, from the mathematical viewpoint, is that it leads to a closed hierarchy of marginal 
equations (or BBGKY hierarchy). We will make use of this opportunity to provide explicit 
computations of the correlations, i.e. of the distance to statistical independence. 

The BCD and CL processes are special examples of a general class of pair-interaction 
processes. The paper will study these processes in the space-homogeneous case. Having 
in mind the special examples of the BDG and CL processes, we assume that the particle 
velocities are two-dimensional vectors of constant norm. However, this assumption could 
be easily waived. The main theorem is that the chaos propagation property is true for 
these pair-interaction processes. The derivation of the BDG and CL kinetic equations 
follow from this theorem. In the case of the BDG operator, we recover the collision 
operator or [2] . The proof of the theorem generalizes some of the combinatorial arguments 
of Mark Kac [26]. 

This result can be seen as paradoxical at first sight. Indeed, the BDG and CL processes 
build-up correlations, in the sense that particles tend to eventually become close to each- 
other (in velocity space). This correlation build-up must be present in any model that 
is to display "swarming" or "flocking" behavior, and it should lead to a breakdown of 
the statistical independence of the particles. The resolution of this apparent paradox is 
seen through an inspection of the time scales. Indeed, the theorem is only valid on any 
finite time interval at the kinetic scale. At this scale, the average number of collisions 
undergone by each particle on a finite time interval is finite and uniformly bounded. If 
is large and one selects a typical pair of particles and a finite time t, none of the particles 
that have collided with the first particles by that time have also collided with any particle 
that has collided with the other particle at time t. If one fixes N, and lets t become large, 
pair correlations do develop. However, if one fixes t and lets become large, the amount 
of pair correlation built up by time t goes to zero as A^ increases. 

At equilibrium, the solution of the master equation converges to the so-called invariant 
density. Due to the special properties of the CL process, it is possible to provide closed 
expressions for the marginals of the invariant density. When the noise distribution is 
appropriately scaled as A^ — ?■ oo, we show that the invariant density cannot be chaotic. 
Indeed, while the single particle marginal density is uniform, the two-particle marginal 
density is not. Therefore, the two-particle marginal density is not a tensor product of 
two copies of the single-particle marginal density, as it should if it would be chaotic. 
Numerical simulations, reported in a forthcoming work [IT], confirm this result. Again, 
it may seem paradoxical that the invariant measure has strong pair correlations and yet 
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chaos is propagated, and again, it is consideration of the time scales that resolves the 
paradox: The time for the particle system to approach equilibrium grows so rapidly 
with that the stationary measure in not relevant to time evolution on the kinetic scale. 

The message of the paper is that the chaos property may be true even for processes that 
seemingly build-up correlations. However, in this case, the correlation build-up capacities 
of the processes under considerations only manifest themselves at scales which are large 
compared to the kinetic scale. To describe these systems at these large scales, kinetic 
theory is not valid anymore, and alternate theories must be devised. So far, the subject 
is widely open in the literature and constitute a fascinating area of research. 

We conclude this section by a few more bibliographical remarks. In |16], an alternate 
kinetic model for the Vicsek system has been proposed. It consists of a nonlinear Fokker- 
Planck equation. It has been derived from a mean-field limit of the Vicsek system in [6] . 
In biological swarm modeling, most of the authors make use of particle (aka 'Individual- 
Based') models (see e.g. [H [131 EH ESI US]) and sometimes, fluid-like hydrodynamic models 
(see e.g. [T3| |35| B2]). The use of kinetic models is more rare. There is a kinetic version 
of the Cucker-Smale model [12]. This kinetic Cucker-Smale model takes the form of a 
nonlinear non-local Fokker-Planck equation which has been rigorously derived from the 
mean-field limit of the discrete Cucker-Smale model in [5]. Kinetic models have also been 
proposed in the context of fish schools [IT] , bacteria and cell motion [221 EH] and ant-trail 
formation [1]. In most cases, their justification is purely formal. The present paper is the 
first step towards a justification of kinetic models in biological swarm modeling. 

The paper is organized as follows. In section [2l we derive the master equations for 
the two examples of biological swarm models that we will consider, the BDG and the 
CL processes. We also provide the definition of the most general pair interaction driven 
master equation. Section |3] is devoted to the proof of the chaos propagation theorem for 
pair interaction driven master equations and its application to the BDG and CL dynamics. 
This proof is modeled in part on the orginal approach of Kac, but with some differences. 
Section H] investigates the invariant measure for the CL process and shows that, under 
some suitable scaling of the noise distribution, it violates the chaos property. A conclusion 
is drawn in section O 

2 Motivation: biological swarm models 
2.1 Two examples 

We consider a population of A^ agents, which to be concrete, we take to be fish in a 
shallow pond, each swimming at unit speed. In this work, we are only concerned with 
the evolution of the velocities, and we neglect the precise spatial location of the fish; that 
is, we assume that all A^ under consideration are sufficiently close to interact with one 
another. 

The shallow pond is essentially a planar domain, and so the individual velocity vectors 
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belong to the unit circle The state space of the system is therefore the torus 



The state of the swarm, or school, is then specified by giving a vector 

V = {vi, . . . ,vn) . 

All agents, or fish, are considered as part of a local cluster, and all are interacting with 
one another. The evolution of v = {vi, . . . , vn) will be modeled in various ways, all based 
on the following general scheme: there is a steady Poisson stream of jump times, at which 
a pair is selected at random from {!,..., A^}, and then these two fish adjust their 
velocities in some way 

(v,,v,)-^iv[,v'^) . (2.1) 

To complete the specification of the dynamics, we need to give the precise rule for updating 
the velocities in (12. ip . Here are the rules we shall consider. 

(1) BDG dynamics: This rule is designed to lead to a kinetic model first investigated by 
Bertin, Droz and Gregoire [2]. The idea is that the pair of agents adjust their velocities 
cooperatively to achieve the same direction of motion, apart from some noise in their 
adjustments. More precisely, define 

Vi . = , ' \ . 2.2 

Now think of Vij G as unit complex number. Let Wi and Wj be two more unit complex 
numbers, chosen independently at random from a probability distribution g{w)(\w on S^, 
and define 

v[ = WiVij and Vj = WjVij . (2.3) 

Regarding as the unit circle in the complex plane, WjVij simply means the product in 
the complex plane of the random variables Wi and Vij, and likewise for the other term. 
Thus, if we write 

Wi = e^® and Vij = then WiVi^j = e'^^+''^ 

so that the noise is additive in the angles. In the case of no noise, Wi and Wj are the 
constant random variable Wi = Wj = 1, and in the case of small noise, they are random 
variables that are strongly peaked around 1. We suppose that g{w) is symmetric; i.e., 

g{w) = g{w*), (2.4) 

where w* denote the complex conjugate of w, and that g{w) is somewhat peaked near 
w = 1. 



5 



Remark 2.1 This rule could be referred to as the "Maxwellian BDG" dynamics, in refer- 
ence to the fact that the selection of the pair {i,j) is independent of their relative velocity, 
like the Maxwellian molecular interaction in rarefied gas dynamics. A more general set- 
ting would make the collision probability of the pair {i,j) depend on their relative velocity 
ViV*, but this will be discarded here, for reasons developed below (see remark WT^) . 

(2) "Choose the leader" (CL) dynamics: In this variant, one of the two agents in 
the pair decides to adopt the other agent's velocity, though it does not get this velocity 
exactly right: The new velocity it adopts is the velocity of the other agent up to some 
noise term. The process is written as follows: if the pair selected is (1,2), and agent 1 
decides to adopt the velocity of agent 2, then the velocity of agent 1 is updated as follows: 

Vi v[ := Zv2 , 

where Z is an independent random variable with values in and probability g. As before, 
we regard as the unit circle in the complex plane, and Zvi simply means the product 
in the complex plane of the random variables Z and Vi. Here again, we assume that g is 
symmetric and satisfies (12 .4^ for the sake of simplicity. We also have in mind that g is 
peaked around 1. 

We use a fair coin toss, modeled by a Bernoulli variable -B, to decide which agent 
adopts the velocity of the other, and then we have the following description of the jump 
when pair (^, j) The selected pair of velocities is updated according to 

v[ = Bv, + {l-B)Zv, 

v] = BZv, + {l-B)vj ^ 

and all other velocities are unchanged. 

There are many other variants on these basic examples, but for now, let us focus on 
these two and seek a passage from this description of the interactions of individual agents 
to an evolution equation for the statistical distribution of the velocities in the system. 
For this we shall employ methods of kinetic theory that have been developed for a similar 
problem concerning colliding molecules in a gas. We shall use a probabilistic approach of 
Marc Kac [26], using a so-called Master equation. 

2.2 Master Equations 
2.2.1 General framework 

We now derive master equations describing the evolution of the probability density for the 
state of the system as it undergoes our stochastic processes. An advantage with starting 
at the microscopic level, i.e., the level of individual agents, is that the modeling is much 
clearer before any large limits are taken. 
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We note that both the BDG or the CL dynamics are Markovian. In general, we 
consider a Markov process on T^v and we denote by G Tat its state just after the kth 
jump. We define its Markov transition operator Q as usual by 

Q^{v) = E{ip{Vk+i) \ V, = v} , 

for any continuous test function (p on T^v- 

Now let Fk{v) denote the probability density of Vk (with respect to the uniform measure 
on Ttv). Then by definition, one has 

JTn 

On the other hand, by standard properties of the conditional expectation. 
That is, 

Fk+i = Q*Fk, (2.1) 

where Q* is the adjoint of Q in L'^(Tn, d^v). 

The next step is to construct a time continuous process which will lead to a time- 
continuous master equation. The state of the process is now a function of time v{t) G T^v 
and the probability density F{v, t) is a function of the continuous time parameter t instead 
of the discrete jump index k. We assume that F{v,t + dt) only depends on F{v,t) and 
not on the past values F{v, s) for s < t. In this way, we can construct a time-continuous 
Markov process. Thus, we assume that between time t and t + dt the probability that a 
collision occurs is 

XQ*F{v,t) dt + o{dt), 

where A is a constant. We assume that the probability that multiple collisions occur 
in the time interval [t, t + dt] is negligible. Thus, since there are particles colliding 
independently, we have, for dt: 

F{v, t + dt)= N\Q*F{v, t) dt+{l- N\F{v, t)) dt + o{dt). 

Equivalently, we can write 

F{v, t + dt)- F{v, t) = NX[Q*F{v, t) dt - F{v, t)] dt + o{dt), 

which leads to the Master Equation in the limit — )■ 0: 

j^F{v,t) = XN[Q*-I]F{v,t), 

where / is the identity operator. This equation must be complemented by the initial 
condition 

F{v,0) = Foiv) 
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where Fq is the initial probabihty distribution. In the remainder, we scale time in such a 
way that A = 1 and we define 

L = N[Q-I] and L* = N[Q* - 1]. (2.2) 

We summarize the previous discussion in the following definition: 

Definition 2.1 The time- continuous master equation associated to a discrete Markov 
process of transition operator Q is written 

j^F{v,t) = L*F{v,t), (2.3) 
Fiv,0) = Foiv). (2.4) 

where L* = N[Q* — I] with Q* the adjoint ofQ and Fq is the initial probability distribution. 

In the next sections, we determine the master equations of the BDG and CL processes 
successively. 

2.2.2 The BDG dynamics 

We state the following: 

Proposition 2.2 The master equation of the BDG dynamics is written Ii2.3\) with 

^* = ^(^2) (2-5) 
and the binary interaction operator Q*^. given by: 

Qh)Fi^) = / ^(^1' ^iv) c/ (viylj) 9 {i^jVlj) ^Vi^^Vj ■ (2-6) 

JT2 



Proof: iFiom. the definition (12. 3 p of the BDG dynamics, we get 

= Elp{vi, . . . , WiVij, . . . , WjVij, ...,vn) 

y2 / ¥^{^1, ■ ■ ■ ,WiVij, . . . ,WjVij, . . . ,vn) g{wi) g{wj) dwidwj 



N{N - 1) f^^ 
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To compute the adjoint of Q, we note that for any probabihty density F on T 



N, 



df 1 . . .dv 



X 



ip{vi, WiVij, WjVij, ...,vn) g{wi) g{wj) dwidwj 

F{vi, ...,v„...,Vj,...,VN)9 [Vivlj) 9 {yjVij) dvidv 

ip{vi, ...,yi,...,yj,...,VN)dvi... dy^, ...dyj..., dv^ 
where we have introduced the variables 

yi = WiVij and yj = wjVij . 
Now changing the names of variables, we finally have 



N 



Q*F{vi,...,vn) 
2 



1) Jr. 



where 



Therefore, we can write 



Vt + Vj 



N 



i<j 



with Ql^ .^ defined by (ES]). Then, using ([22]), we get eq. ( 1231) . 



Remark 2.2 It is useful to note that the adjoint L of L* which corresponds to the Markov 
transition operator, is defined by 



(2.8) 



with 



Q{i,j)^i^) = / V^K' • • • ) uJi^ij, • • • , WjVi,j, • • • , vn) giw-i) g{wj) dwidwj 
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2.2.3 The CL dynamics 



We now derive tlie master equation for the CL model. We introduce the notation 
(f 1, . . . ,Vi, . . . , Vn) for the n — 1 tuple formed by removing Vi from v and 



[F]-^{vi, . . . ,Vi, . . . ,vn) ■■= / F{vi, . . . ,VN)dvi 
for the marginal of F obtained by integrating Vi out. We show the 

Proposition 2.3 The master equation of the CL dynamics is written Ii2.3\) with L* given 
by 1^2. 5\} and the binary interaction operator Q*^- by: 

Qh)^i^) = I Imvi, ...,Vi,...,VN) + [F]j{v,, ...,v,,..., Vn)] g{v*vj) . (2.9) 



Proof: From the definition (12.51) of the CL process, the Markov transition operator Q is 
given by: 



Q<^{v) 



,ZVj,...,Vj,...,VN) 



+ip{vi, . . . , Wi, . . . , zvi, . . . , Vn)] g{z) <lz . (2.10) 
To compute the adjoint of Q, we note that for any probability density F on T^, 



F{vi, ...,vn) 



^{vi, zvj, ...,Vj,..., VN)g{z)dz 



dvi . . .dv 



N 



[F]-{vi, . . .,Vi,. . .,vn) / ^{vu • • • , zvj, . . . ,Vj,. . . ,vn) g{z) 

_i JSi 

d^; dwi . . . df j, . . . , df TV (2-11) 

We next introduce a new variable 

yi = zvj (or equivalently z = v*yi ). 

Evidently dz = dy.i. Additionally, since Vi has been integrated out, it has disappeared 
from (12.111) . Therefore, we can change the name y^ into Vi without any confusion. We can 
then rewrite (12.111) as 



F(vi, ...,vn) 



ip{vi, . . . ,zvj,. . . ,Vj,. . . , VN)g{z)dz 



dvi . . .dv 



N 



[PhiVi, ...,Vi,...,VN) g{v*Vi) Lp{Vi, ...,Vi,...,Vj,...,VN) 

dvi . . . dvi, . . . , dvN ■ (2.12) 

Using this formula and its analog for i and j exchanged, we see that the master equation 
for the CL dynamics is given by (12.21) . (12.51) with Q{ij) given by (12.91) . ■ 



10 



Remark 2.3 Again, we note that the adjoint L of L* is defined by i\2. 8\) with Q{ij), the 
adjoint ofQ^-y given by: 



Q(ij)Vi.v) = / (fivu . . . , zvj, . . . , uj, . . . , VN)giz)dz . 

2.3 Extension: Pair-Interaction driven Master Equation 

The master equations of the BDG and CL dynamics are two examples of a class of master 
equations which we will call 'Pair Interaction driven Master Equations', defined below. 

Definition 2.4 (Pair Interaction Driven Master Equation) A pair interaction driven 
Master equation is an equation of the form 

^F{v,t) = L*F{v,t), 

describing the evolution of probability densities on some product space with elements 
V = (t>i, . . . , vn) where 

l* = nY,p,,{v){qi^^^-i) . 

The operators Q{ij) are Markov operators on functions on Xn such that Q(^ij)ip = ip 
whenever ip does not depend on either Vi or vj. The pair selection probabilities Pi.j{v) are 
such that Pi J (v) > and 

^Pijiv) = 1 . 

i<j 



We have given two examples already: the CL and BDG master equations. The Kac 
Master equation [26] is another example and is described below. 

Example 2.1 The Kac Master equation. In this example, X^r is the sphere in of 
radius VN, 



J — TT 

RijfiV = {vi,V2, . . . , COS 9vi + sin 9vj, . . . , — sin 9vi + cos 9vj,..., Vn), 
2 

~ N{N - 1) 

and p is a probability density on S^. The operators Q(ij) in the Kac model are self adjoint 
with respect to the uniform probability measure on the sphere S^~^, which is therefore the 
invariant measure for this process. In other words the Kac process is reversible meaning 
that is satisfies detailed balance.- if you saw a movie of the process running backwards, 
there would be no clue that it was running backwards. 
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By contrast to the Kac model, the BDG and CL models do not have detailed balance 
property and time reversibility. If you ran the movie backwards, you would see pairs of fish 
with similar velocities changing them to differ in a random way. For these non-reversible 
processes, it is not so easy to determine the invariant measure, though it will exist and 
be unique for each for our processes under mild assumptions on the noise distribution 
g. In section m it will be possible to determine the marginals of the invariant density of 
the CL dynamics in closed form. However, this simplification is not possible for the BDG 
dynamics. 

In the next section, we show that pair interaction driven master equations with uniform 
selection probabilities Pij = 2/{N{N — 1)) do have the propagation of chaos property, 
and therefore, satisfy a kinetic equation at the kinetic time scale. However, in section HI 
we show that the equilibrium density of the CL dynamics cannot satisfy the propagation 
of chaos property, meaning that this property may break down at larger time scales. 

3 Propagation of Chaos 
3.1 Definition 

To pass to a kinetic description, and then on to a hydrodynamic description, the key step 
is a propagation of chaos result. That may seem unlikely in the cases of the CL and 
BDG dynamics which are expected to build pair correlations. However, the time scales 
at which pair correlations built up may be longer than the kinetic time scale at which 
a kinetic model is expected to be valid. In the present section, we shall see that chaos 
is propagated in both the BDG and CL models at the kinetic time scale. In section HJ 
we prove that the invariant measure of the CL dynamics is not chaotic; it exhibits pair 
correlation. These two observations are not self-contradictory since propagation of chaos 
holds only on a finite time scale while the invariant measure is reached as time tends to 
infinity. 

Definition 3.1 (Chaos) Let X^r be the N-fold cartesian product of a polish space X 
equipped with some reference measure fi. Let f be a given probability density on X . For 
each N G N, let F(7v) be a probability density on Xn with respect to /i*^^. The sequence 
{F[n)}n(^n of probability densities on Xm is /-chaotic in case 

(1.) Each Fi^i^^ is a symmetric function of {vi, V2, . . . , v^} 

(2.) For each fixed k, and any bounded measurable function on M^, 



Kac [26] proved that the semigroup e associated to Kac's master equation propagates 

chaos. More precisely, Kac's Theorem is stated as follows: 




k 
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Theorem 3.2 (Propagation of chaos) Let {-F(jv)}AfGN be f -chaotic. Then the family 
|gtA^(Q-^)^^^^ j^^pj f{-,t) -chaotic where f{v,t) is the solution of 

^{v,t) = Q{fJ){v,t) with f{v,0) = f{v), (3.1) 



with 



Q(/, /) = 2 / r [/(t;')/K) - f{v)fH] p{e) dedw 



and 

v' = cos 9v + sin 6w, w' = — sin 6v + cos ^ . 

Eq. (13. ip is called the Kac-Boltzmann equation. In this section we prove a propagation 
of chaos result valid in the general class of pair interaction driven Master equations. We 
shall use this result to discuss the kinetic limits of the BDG and CL dynamics. 

3.2 Propagation of chaos for pair interaction driven master equa- 
tions 



Consider a general Master equation 



= L*F, (3.2) 



for a probability density F on T^v of the form 

L*F = N{Q* - I)F = J2(Qh) - ^)^' (3-3) 

i<j 

where Qij is a Markovian operator acting on F through Vi and Vj alone. The goal of this 
section is to prove the following: 

Theorem 3.3 Let {FQ^}N^fq be f -chaotic. Then for each t > 0, the family of marginals 

{e*-^*Fg^''}ArgN associated to eq. \3. where L* is a pair-interaction operator of the form 
( 15*. 5]) . is f{-,t)- chaotic where f{-,t) satisfies the following Boltzmann equation: 



^^f{v, t) = 2 Ql^,,)f^'{v, w) dw - f{v, t) 
associated to the initial condition f{v,0) = f. 



(3.4) 



Before proving Theorem 13.31 we make some preliminary comments. Let initial data Fq^n 
be given, and let us compute the evolution of -F^^^ the single particle marginal at time t. 
For any test function ip{v) of the single coordinate v e we have 

[ ip{v) F,^'\v) dv = I ^{vi)e'''''Fo^N{v)dvi...dvN 

e*^(p{vi) Fo^n{v) dvi... dvN ■ 

V 
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A similar relation holds for the two particle marginal and so on. So, to study the evolution 
of low dimensional marginals, it is helpful to understand the behavior of expressions of 
the form e^^ipiv) when ip{v) depends on only finitely many coordinates in v. 
It is clear that in general, for a bounded continuous function ip on T^v, 

IIL^IU <2Ar||^||^ . 

However, if ip depends only on t^i, . . . , w^, tighter bounds are valid. This is because 

hi > k =^ Qii,j)¥^ = , 

and so 

2 2 ^ ^ 

i<j i=l j=i+l 

and thus, as soon as A; > 1: 

||^<^||oc < k{N - ^^) 2||^|U < 4A;||^|U • (3.6) 

We can now state the following fundamental lemma: 

Lemma 3.4 Let ip he a function depending only on vi, . . . ,Vk- We can regard ip as a 
function on Xjsi for each iV G N, N > k. Then, the power series 

^'V^Ef^V (3.7) 

k=0 

converges absolutely in L°° , uniformly in N ^N* and t G [0, T] for any T < 1/4. 

Proof: Consider first the case in which ip depends only on one variable. Without loss of 
generality, owing to the permutation symmetry of the problem, we can set this variable 
to Vi. Then from (13.51) . Lip is an average of functions depending on only two velocities. 
Likewise, L'^ip is a combination of terms depending only on three velocities and so on. By 
what we have noted above, we can expect the following formula: 

Ili^Vll <4'A;!||v.||oo . (3.8) 
To show (13. 8p we prove that L'^ip is of the form 

where the set Sk is a set of multi-indices s = (1, Si, . . . , s^), such that 

k 

Card Sk < k\Y[{N -^-^). (3.10) 

i=i 
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The function if)^^'^ depends only on the k + 1 variables (wi, Wg^, . . . , Vg^) and satisfies 

llV^f ||oo<2'l|y^|U. (3.11^ 
Of course, (13. 8 p results from f l3.9p . f lS.lOp and (13. lip and the remark that 



1 ^ 



iN 



We note that, if A; + 1 > A^, some of the variables (fi, f^^, . . . , Vg^.) may be identical, but 
the arguments below are still valid in this case. 

The proof of ( 13. 9p is by induction. For /c = 1, using ( 13. Sp . we have 



2 ^ 



i=2 



Therefore, letting 5*1 = {2, . . . , A^} and ^/'i^'' = ip^p = {Q{ij) — I)^, we can write Lip 
according to formula ( 13. 9p . Clearly, Card S'l = — 1 in accordance to ( 13.10p . Finally, 
by the Markov property, we have ||'?/^]^''||oo < 2||</)||oo, which is consistent with (I3.1ip . 
Therefore, (13. 8p is proved for k = 1. 

Now, we assume that (13. 8p is true for k and try to deduce it for k + 1. By the induction 
hypothesis, ijjf^ depends only on A; + 1 variables so ( 13. 5 p applies and we have 



^ seSfe 

^ ^ seSfe m£s,e>m 



The expression m & s means that m takes any of the indices {1, Si, . . . , Sfc} present 
in the multi-index s. We know from the computation in formula (13.60 that there are 
[k + 1)(7V — such pairs (m,£) for a given s. Defining S'fc+i as the set of the so- 
constructed multi-indices {s' = {s,i)}, we can write 

with 
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Clearly, ipl'^'^^^ is a function of the k + 2 variables (wi, w^^, . . . , Vsf.,Vm) and we have 

k + 2. 



Card^fc+i = {k + l){N 
= {k + l){N 



2 

k + 2 



Card Sk 

k 

k\ W{N 



J + 1 



fe+i 

{k+l)\\{{N 



Finally, we have, using the Markov property 

||^?+'^||oo<2||^f |U<2^+1||V^|U. 



(3.13) 



(3.14) 



Now, collecting fl3.12p . fl3.13p . fl3.14p shows that the induction hypothesis is valid at rank 
k + 1. This proves (13.81) . 

Using (13. 8p . we deduce that 



k\ 



<m'M 



uniformly in A^, and so for t < 1/4, we have the uniform absolute convergence of the 
series (13. 7p . 

If now (y9 is a function p variables, p > 2, formula (13. 8 p is changed into 



p + k — 1 
k 



< A'^kl 



{p-i)\ 



< Cp4}'k\{k + lY~^\\^\ 



(3.15) 



(3.16) 



The proof follows the same lines as above. The only thing to note is that Card is now 
changed into 

k . 

Card^fc = n(P + J-l)(A^-^), 
i=i 

the other expressions remaining identical. Then, we get 



k\ 



<Cp{u)\k + iY-^M 



(3.17) 



uniformly in A^. For any given p >2, the right-hand side of (I3.17P is still the general term 
of a convergent series for t < 1/4. Therefore, the series (13. 7p is still absolutely uniformly 
convergent for t < 1/4, which ends the proof of Lemma [3.41 ■ 
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Remark 3.1 iFrom the last proof, we note that the series (3.1) is not uniformly con- 
vergent with respect to p. 

We now have the following 

Lemma 3.5 Suppose that Fq^n is a symmetric probability density on T^, and suppose 
that ip^^^ depends only on vi, . . . ,Vk and is L°°. Define ip^^'^^^ by 

k 

<^('=+i)K, . . . , Vk+i) = 2 5^(g(.,fc+i) - /)<^('n^i, • • • , ^fc) • (3.18) 

i=l 

Then, if N > k + 1, we have 

[ Fo,NLip^''^dv,...dvN= [ Fo,Ni^^''^'^ + ^^''^'^)dvi...dvN, 
where (p^'^^^^ only depends on fi, . . . ^Vk+i and is such that 



(3.19) 



U''^''\U<Q^-^^\W'''\U. (3.20) 



Proof: We compute, using (13. 5p and that N > k + 1: 
I Fo,jv Lip^^^ dvi... dvN = 

k N ^ 

Fo,N{Qi^J)-I)^^'^dv,...dvr, 

i=l j=i+lJ'^N 



2 

w 



k 



i=l -^"^N 

+ E / Fo,NiQi^,)-I)p>^'^dv,...dv 

A ^ A^l^ T A/- 



where we have used the symmetry in the second equality, and where 



i=l 



N _ 

i<j<k 
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This shows fl3.19p . Now, from the Markov property of Q(ij) and from fl3.2ip . we get fl3.20p . 
which ends the proof of the Lemma. ■ 



Now consider any yj^*"^ depending only on Vi, . . . ,Vm- Since if Fq^n is symmetric, so is 
each {L*)^FQ^i\f. Therefore, we can repeatedly apply the previous lemma and so on. Using 
IKTE^ . we inductively define . . . (^("^+*^). We note that, if does not 

depend on A^, neither does (p^"^^''\ Therefore, the functions tp^"^^''^ are good candidates 
to express what happens in the limit — > oo. However, for a given N, formula (13.191) 
is only valid until m + k < N. A special treatment is required for indices k such that 
m + k > N. The contribution of the remainder cp^'''^^^ needs also to be estimated. These 
are the goals of the following Lemma. 

Lemma 3.6 We assume that {-Fo.tv} is f -chaotic (see definition \3 . 1\) . Then, fort < 1/4, 
we have: 

hm / Fo TV e*^</?('") d^;i . . . dw^v = 

oo ^ „ /k+m \ 

= E / n ^(^'i) ^^'"^'H^i, • • • , ^fe+rn) d^i . . . di;fe+„ . (3.22) 

fc=0 ■ -^"^k + m. \j = l J 



Proof: j'From Lemma [3.41 we have 



with 



lim / Fo^N e*^</?('") dvi... dvN = lim lim S^{t) 



^nit) = "E 77 / ^0.^ ^V^'"^ dvi... dvN , (3.23) 

Indeed, according to Lemma 13. 4[ the series (I3.23P converges absolutely, uniformly with 
respect to N and we can interchange the if — )■ oo and N oo limits. Now, for k+m < N, 
according to the inductive definition of (^('"+'^) and to (I3.19p . we can write: 

/ Fo,N ^V^™^ dvi... dvN = [ Fo,iv f v^('"+'=) + V L*^-^(^("^+^M dVi... dVN. 
Taking an index H such that H + m < N, we have: 

S^it) = E 77 / ^0.^ '^^"^'^ • • • + / ^0'^ ^^'^^"^ • • • (3.24) 



k=0 

with 



k=0 ' j=l 



Now, taking successively — oo then i7 — )■ oo, we show that the first term of (13.241) 
tends to the right-hand side of f l3.22p and the second one tends to zero. 

We start with the first term. By the chaos assumption, we have, as — )■ oo: 



hmV-/ Fo,;v</^^'"+'Mt;i...dt;^ 



fc = ■ "^^m+fe \i = l / 

We show that the series at the right-hand side of (13.251) is absolutely convergent, uniformly 
with respect to t in any interval [0,T] with T < 1/4. The proof follows the same lines 
as that of Lemma 13.41 and we only sketch it. Using (I3.18P and the Markov property of 
we have 



^(m+fc-l) 

m + k — 1 
k 

< Cm^'' k\ {k + ir-' y^""^ 



< rk\ 



II 

V II CO 



(3.26) 



Therefore, since n^=r fi^j) a probability: 



//k+m \ 
-'m + fe \j=l / 



<cu'^t)'ik + ir-'y^"'^\U 



This is the general term of a convergent series which converges uniformly as stated above. 
We deduce that 



lim lim V - / FoN^^""^''^dvi...dvN 

H^oo N^oo ^ k\ T 
k=0 



t 

'hi 



k=0 



^ k+m 

n /(^^-^ 

1+fe \jr = l 



N: 



which is the right-hand side of (I3.22p . 

We now consider the second term of ( 13.24^ . First, using (13.15^ with the pair (fc,p) 
replaced by {k — j,m + j), we get 



,,^k-j~(m+j)u ^Ak-j i^ + k ^V- ||^(™+i)|| 
\\J^ y-> oo ^ ^ / , . 1 M IIV' lloo 

(m + J — 1)! 

Then, (I3.20p with k replaced by m + j — 1 yields 



(3.27) 



(3.28) 



19 



Finally, using f l3.26p with A; = j — 1, we obtain 



^('"+^-^)||oo<4^'-^%^^^-r#||^('")||oo. (3.29) 
[m — 1)1 



Collecting (K27} . (Km and (Km leads to 



and: 



" ^ iioo _ 2 (m-1)! N-1 ' 



^ k-j~im.+j). < 3 {m + k-l)\ jm + k) 
^ ^ Hoc S 2^ 2(iV-l) ' 



which finally gives: 



m+l 



i=i 

Therefore, we have 



(m)| 



I OO) 



/ Fo,N R^^^""^ dv,... dvN <^Crn ^^""^ ||oo ( V(4t)'^ (A; + 1)"+^ 

>^Tiv Vfc=0 

and deduce that 

lim lim / Rh^^"^^ dwi . . . dv^ = 0, 

which ends the proof. ■ 
We now have all the elements to prove Theorem 13.31 

Proof of Theorem 13. 3t We need to show the existence of a function of a single velocity 
variable f(^v,t) such that 



lim / Fo,^eV"'^clt;i...dt;;v= / \flf{vj,t) 



dvi... dvN, (3.30) 



for all functions (p^"^^ of m velocity variables (wi, . . . , Vm) and all m < 1. 

First applying this definition with m = 1 defines f{v,t) (if it exists), by duality: 

/ f{v,t) ip{v) dv = \im / FoN e^^ipdvi...dvN- (3.31) 



20 



Now, applying f l3.30p with m = 2, </3^^^(i'i, 1^2) = ''7''^''(^'i)'^''"'^''('^2), and using (13.311) leads to 
lim / FoNe^^if'-'^^ dvi . . .dvN = 

lim / Fo TV e*^?]^^) dvi . . . dtJ^v ) ( lim / Fq n e*^^^^^ dvi . . . dvN ) ■ (3.32) 



Reciprocally, (13.321) implies (I3.3QP for m = 2 and general functions ip^'^\vi,V2) by the 
density of linear combinations of tensor products. Finally, the same argument applied to 
arbitrary m would require to check the property for (p^"^^ equal to the tensor product of 
m one- dimensional functions. This proof follows closely the proof for the case m = 2 and 
will be omitted. So, we now focus to the case m = 2. 

We first consider the case of small t < 1/4. Then, using (I3.22p . eq. (I3.32p is equiva- 
lently written: 



/ n/(-.)K^'''-'^d.,...d.^^2= 

K=0 ' \j = l J 

nM)h^'"^^d.,...d..,, 



X 



X 



E^ / nM)h^'"'^d.,...d.,^, 



We start with the right-hand side of this formula, denoted by Tt. By distributing the 
various terms in the product, we can write 

^-t%i (nV)) (£(::; 

r/("^+i)(i;i, t;3, . . . , W2)e^''-")+'(^2, W3, • • • , ^^^+2) j d^;i . . . dvK^-,. 

Therefore, the result is proved if we show that for any symmetric function of F(wi,U2, 
• • • ,'^^-+2), we have 

^ f K \ 

• • • 5 Vk+2) 



[ E ( f )^^™^^H-i,-3,...,w2)e(^-")+^(^2, 

\m=0 ^ ^ 

F{vi, V2,..., VK+2) dVi... dVK+2 = 

= / ip^^^^\vi,V2, . . . ,Vk+2) F{vi,V2, . . . ,VK+2)dVi . . .dvK+2 ■ (3.33) 



21 



Indeed, at the initial step K = 1, by direct computation we get: 

(p^^\vi, V2, Vs) F{vi,V2, Vs) dvi ...dV3 = 

T3 

Then, f l3.33p is easily proved by induction, using elementary properties of the binomial 
coefficients. This shows that e^^* Fq^n is /(■, t)-chaotic for small t < 1/4. 

It remains to show that f{v,t) is a solution of ( 13. 4p . Again, assuming small t < 1/4 
and Using f l3.22p with m = 1 and f l3.3ip . we get 



nv,ty'Hv)dv = j:j^ n/(' 

"^^1 k=0 -^^k+i 



j) </3^''^^n^'i> • • • 5 ^fc+i) dvi... dvk+i . 



The convergence of the series at the right-hand side is uniform for t < 1/4. So, we can 
differentiate this formula with respect to t and obtain: 

^{v,t)^^'\v)dv = 

Ti <^ 

= E / n /(^^O ^^'^'H^i, • • • , dv^... dv,^2 . (3.34) 

Applying (13.22P and ( ]3.30p with m = 2, we can re- write the right-hand side of ( 13.34p and 
get 



l^^(v,t)v<^'\v)dv = 1^ (^J{v„t)^ ^^'^dv,dV2 



2 / fivu t)f{v2, t) (g(i,2) - IW^'^ dVi dV2 

{Qli,2)-m^f)¥>^'^dv,dv2 , 

which is the weak form of (13. 4p . This shows that f{v,t) is a weak solution of (13. 4p for 
small t < 1/4. 

In the last step of the proof, we need to remove the restriction on t < 1/4. However, 
since this bound is independent of the initial data, we can partition any interval [0, T] by 
intervals [tfc, tfc+i] of length tfc+i —tk < 1/4 and apply the result on each of these intervals 
with initial data Fjsfitk) and f{v,ti). This shows that e^^* Fq^n is /(■, t)-chaotic on any 
finite-size interval [0,T] where f{v,t) is the solution of (13.40 . This ends the proof. ■ 
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Pi,j = 1^771^7 —, (3.35) 



Remark 3.2 In Theorem \3.3\ the assumption that the pair selection probabilities Pij{v 
are uniform, i.e. 

2 

N{N- 1) 

is crucial. Indeed, we need at least that 

(i) there exists a uniform constant C (independent of N) such that 

C 

Pi,j < 



N{N-1) 



This is required to get the fundamental property ( IJ. 6]) . with the constant Ak replaced by 
ACk. 

(a) pij = pij{vi, Vj). In this way, the key fact in the proof of Lemma \3.4\ that Lip, for ip 
depending on only one velocity, is an average of terms only depending on two velocities 
and so on for L^ip remains true. 

(Hi) Pij{v,w) = p{v,w) is independent of to preserve the permutation symmetry of 
the problem. 

We see that these three properties together imply that 35\) . Indeed, maintaining that 

5^p„ = l, (3.36) 

i<j 

for all velocity configurations necessitates that p is a constant, and the normalization 
constraint Ii3.36\) leads to Ii3. 35\) . 

Remark 3.3 Theorem \3.3\ can be extended to interaction processes involving multiple 
interactions as soon as the number of particles involved in an elementary interaction is 
finite and bounded independently of N . For instance, it will hold with a ternary interaction 
process, in which any interaction involves triples of particles. More generally, it will hold 
with a p-fold interaction process where any interaction involves exactly p particles. As 
long as the interactions involve a finite number p of interactions, with p constant or 
bounded by a constant P independent of N, the combinatorial arguments which have been 
developed above can be extended. Again, the master equation must combine the elementary 
interaction operators by means of uniform selection probabilities for the same arguments 
as those developed in Remark \3.2[ Boltzmann operators with multiple interactions have 
been previously considered in 

Remark 3.4 74s a by-product of Theorem \3.3\ we have the existence of a weak solution 
of the nonlinear problem ^3.4^ . Additionally, the proof is constructive. However, the 
Theorem does not prove uniqueness, as another solution could exist not being related to a 
solution of the corresponding master equation. 
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3.3 Application to the BDG and CL dynamics 
3.3.1 The BDG dynamics 

In the BDG case, thanks to (12.61) . we have 



and thus 



To 



Thus defining 

g+(/, /)(^) = / /(yi) /(z/2) 2) dyi dy2 , 

the one particle marginal /(■, t) such that {e^'^'-fo at } is /(■, t)-chaotic satisfies the kinetic- 
type equation 

^{v,t) = f{yi)f{y2)g{vyl2)dyidy2-fiv,t) , (3.37) 
with/(i;,0) = /(i;). 

3.3.2 The CL dynamics 

In the CL case, thanks to (12. 9p we have: 

1 



QK2)f^\v,w) = -{f{v) + f{w))9{v*w) 



and thus 



Qii,2)f^\^,^)dw = -i^fiv) + jj{w)g{v*w)dwj=-[f{v) + f^g{v)] 
where -k denote the convolution. Thus defining 

Q+UJ){v) = \[f{v) + f^g{v)] , 

the one particle marginal f{-,t) such that {c^'^'Fo^at} is /(■, t)-chaotic satisfies the kinetic- 
type equation 

^{v,t){v,t) = ^[g^f(v,t)-f{v,t)], (3.38) 

with f{v,0) = f{v). 

In this treatment, we have assumed that g is independent of A^. But if we let the 
variance of gN go to zero with A^, we find 

lim gN-*^f{v) = f{v) , 
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and then we have 

|(M)^0. 

That is, chaos is propagated, but nothing at all happens on the kinetic time scale. On 
a much longer time scale, correlations develop and a new approach is needed to describe 
the bulk limit. This is what is shown in the next section by investigating the invariant 
densities, i.e. the equilibria of the master equation. 



4 The Invariant densities F^^^ for the CL dynamics 
4.1 Preliminaries 

Both the BDG and CL processes are clearly ergodic as long as g is continuous, say, and 
so for each there will be a unique invariant density F^o, i.e., a unique density F^o with 
Q*Foo = Foo- Since the process is symmetric under permutations of the variables, it is 
clear that Foo will be symmetric. It is not easy to write Foo down in closed form. However, 
in the case of the CL dynamics, a very special property is true, namely the hierarchy of 
equations for the marginals (or BBGKY hierarchy) is closed at any order. This special 
feature will provide more information on the marginals of Foo. In particular, we show 
that under some specific scaling of the noise with respect to A^, the invariant measure is 
not chaotic. This is not in contradiction to Theorem 13. 3[ since it is valid for fixed noise 
and on finite time intervals. 

The invariant density for the CL dynamics is the function F^o which cancels the right- 
hand side of (12.51) with Q*^ij) given by (12.91) . Therefore, it satisfies 



+ [Foo]^{vi, ...,Vj,...,VN) g{v*Vj) . (4.1) 



While it is not easy to write F^o down in closed form, we can at least say what Foo is not 
In general F{v) = 1 does not solve Q*F = F, i.e., Foo is not the uniform density. Indeed, 
if we replace Foo by 1 on the right hand side of (14.11) , we find 



N(N 



r^(^j^i) 



This will equal 1 for all v if and only if g{z) = 1 for all z. However, it is easy to see that 
for fixed smooth g, and large N, 



^iV(iv3T)^^^'.-^^-^' 

i<j ^ ^ 

with high probability if v is selected at random, uniformly on Tat, and so one can expect 
that for fixed g, the invariant density Foo becomes more and more uniform as increases. 



25 



However, a non-uniform invariant density can be found in the large limit if one 
includes some dependence in the noise density g in such a way that it more and more 
closely approximates a 6 function. This is perhaps justified in the context of biological 
modeling: if the population is small, one fish may be less interested in carefully mimicking 
his neighbor than when the population is large. One can imagine that the larger the group, 
the more important it is to follow behavioral rules closely. However, biological data are 
needed to support this claim. 

While it does not seem easy to write F down in closed form, it is possible to obtain 
analytical expressions of its marginals. This is the aim of the next section. 

4.2 Marginals 

For any symmetric density F and any m = 1,2, — 1, define the m- variable marginal 
density F*^"*) on by 

F^'^\vi,. . . ,Vm) = j F{vi,...,Vm,Vrr,+l,...,VN)([v^+i...dVN- 

For typical binary collision process master equations, the evolution of the m-variable 
marginal density depends on the (m + l)-variable marginal density, and one gets a hi- 
erarchy of evolution equations, the so-called BBGKY hierarchy. To break the hierarchy 
at some manageable level (e.g. m = 1 in the usual kinetic theory case) in the large N 
limit, one generally needs special assumptions on the initial data, such as the /-chaotic 
property, for a suitable /. One typically runs into the same hierarchy problem when 
trying to compute the marginals of the invariant density, except that in some cases, one 
has Q = Q*i and then the invariant density is simply uniform. 

For the CL dynamics, the hierarchy breaks itself. Before stating the result, we recall 
that the Fourier series of a given function (piw), G is defined by 

4){k) = / w~^(f){w)dw, keZ. 

We have: 

Proposition 4.1 Let f!;^^ and Fc^'' be the one and two-variable marginal invariant den- 
sities of the CL process. We have: 

(i) F^ is uniform, i.e. 
(a) f!^ is given by: 

(4.2) 



with J-'{w), w E 'B given by its Fourier series 

1 



N -1 



1 9{k) 
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or equivalently 



7 



oo 



N -2 



N-2 
N-1 



(4.3) 



Remark 4.1 Since Yl'eLi {w^Y = N — 2, Eq. ( [^.^ defines T as an average of convo- 
lution powers of g . 

Proof: We begin with F^-*. First of all, it is easy to see that for alH > 1 and all j , then 
[Foo]i(wi, ...,Vi,..., VN)g{v*Vi)dv2 ■ ■ ■ dVN = F^\vi) , 

-1 

Therefore, by integrating (14. ip with respect to (^2, • • • ,vn) and using (12 ■4p . we have: 



I — ttY] [ [F^]^{vi,...)g{vlvj)dv2---dvN 

N-1 



N 



F^\v,). (4.4) 



Then, using the permutation symmetry of [-Foo]t, one finds for each j > 2, 

[i^oo]T(^^i, • • • ) 9{v*iVj) dv2--- dvN = / F^\vj) g{vlvj) dvj = * g(vi) 



where the * denotes convolution on S^. Substituting it into (I4.4p . we find 

F^\v,) = * g{v,) + ^F«(^i) , 

which reduces to 

F!;^\v,) = F^U g{v^) . 
The only solution of this equation for any density g is the uniform density f!^\vi) = 1. 

(2) 

We now turn to Fso {vi, V2). Again, it is easy to see that if i,j > 2, 



Tjv- 



and 



[Fookivi, V2, . . . ,Vi, . . . , vn) g{v*Vi) dw3 ■ ■ ■ dVN = F^\vi, V2) 



[i^oo]7(^^i,^^2, ...,Vj,.. .,VN)g{vlvj)dv3 ■ --dyN = F^\vi,V2) 



Next, for z = 1 and j = 2, using the permutation symmetry of [-Foo]t and (12 ■4p : 
[F^]iivuV2, ...,vn) 9iv*vi) dv-i--- dvN = F^\v2) giv*vi) = g{vlv2) 



Tjv- 
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and 

/ [^oo]2(^'i,^2, • • .,VN)g{vlv2) dv-s ■ --(ivN = Fj;^\vi)g{vlv2) = g{vlv2) . 
Finally, for i = 1 and j > 2, we have 

[Foo]i(tri, W2, . . . , vn) g{v*vi) dv3--- dvN = / F^\v2, Vj) g{v*vi) dvj , 

Tn-2 

and 

/ [Foo]j{vi,V2, . . . ,Vj, . . . ,vn) g{vlvj) dvs ■ ■ - dvN = F^\vi,V2) , 
and for i = 2 and j > 2, 

/ [Fochivi, V2, Vn) g{v*V2) dv^ ■ ■ ■ dvN = / F^\vi, Vj) g{v*V2) dvj , 

and 

/ [^oo]j(t'i, V2,...,Vj,..., Vn) g{vlvj) dv-i ■ --dvN = F^\vi, V2) . 

Tjv-2 

We can now compute the two variable marginals of both sides of (14. ip . and we find 



F^\v,,V2) 



{N-2){N-3) ,2 



N{N-1) N{N-1) 



FB\v,,V2) 



N{N-1) 



g{vlv2) + .."fl, '^,^ H{vi, V2), 



where 



H{v,,V2)= / F!^\v2,z)giz*v^)dz+ / f!^\v„v,) giz*V2) dz . 
This simplifies to 

F^\v„V2) = j;^g{vlv2) + ^^^^(^'1.^2) . 

Fourier transforming both sides, we have 

H{h, k2) = F^\h, k2){g{h) + g{k2)) 

and so 



N -2 f g{k^)+g{k2) \ 



iV- 1 V 



A^- 1 



g{h)Sku- 



k2 



where Sij is the usual Kronecker symbol. It follows that F^\ki,k2) has the form 

J^{ki)Sk^-k2, i.e. that f!^^ has the form f!^\vi,V2) = J^{vlv2) with defined by its 
Fourier transform according to (14 ■2p (owing to the evenness of ^f, a consequence of (I2.4p ). 
This ends the proof of proposition 14.11 ■ 
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4.3 Noise scaling 

The reason for considering a scaling of the noise intensity is the following. For fixed 
g G lim^^oo 9*^ = 1, and since for large A^, most of the weight in the average is on 

large values of i, T will be nearly uniform for large values of A^, and the correlations are 
washed out. Therefore, we recover here that the invariant density is nearly uniform for 
large values of A^, a fact which has already been noticed (see section |4]T]) . 

But if g is taken to depend on N itself, this need not be the case. As a typical example, 
we can consider g(z*w) to be the kernel of e^^^ on S^; i.e., the heat kernel on at time 
1/A^. Then 



9N{k) :-- 



-k-'/N 



(4.5) 



We now state the 



Proposition 4.2 Suppose the scaled noise intensity g^^ is such that 

lim {N-2){gr,{k)-l) := ^{k) 

N^oo 



(4.6) 



exists and is non trivial (i.e. not equal to the Kronecker 6k,o). Then, the corresponding 
correlation J-'m associated to gN through Ili4.2\ ) or ( [^.g| ) satisfies: 



lim TN{k) := J^oo{k) = ttv- 

N^oo 1 — 7(A;) 



Proof: Using (7 at in place of g in (14. 2p . we find 



J'Nik) 



1 



A^- 1 
1 

A^- 1 



QNik) 
9Nik) 



1 - 



A^- 2 



A^- 1 
1 A^-2 



A^- 1 A^- 1 



(1 + (Mk) - 1)) 
{9N{k) - 1) 



-1 



Mk) [l-{N-2){Mk)-l)r' 



The result follows from inserting (14. 6p . 



(4.7) 



Example 4.1 The heat kernel In this case, we have 

l{k) = -k\ ^oo{k) ^ 



1 + k^' 

Hence, the correlation function is a Lorentzian in Fourier space. 
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With this noise scahng, the two-variable marginal invariant density of the A^-particle 
CL process F^^{vi,V2) is such that 

(2) 

-^oo,7v(^'i' ^2) J^oo{vlv2), as N ^ 00, 

where Too is not the uniform distribution. Therefore, non-trivial correlations remain in 
the large limit and in particular, the invariant density {-Foo,Ar} is not chaotic. This 
result is in marked contrast to the case studied by Kac, in which the invariant density is 
the uniform density on the sphere S^^^{\/N). This family is well known to be G-chaotic 
where G{v) denotes the centered unit Gaussian on M. The lack of chaos in the invariant 
density might seem to be a strong obstacle to propagation of chaos. But as we have seen 
in Theorem 13. 3[ this is not the case. As already noticed, there is no contradiction between 
this two seemingly paradoxical results. Theorem 13.31 is valid under fixed noise and on a 
finite time interval. By contrast, the lack of chaos property of the invariant density for 
the CL model is shown under iV-dependent noise intensity and in the infinite time limit. 
But these two results show that the chaos property can be valid for a finite time interval 
at the kinetic time scale and be lost at larger times. 

5 Conclusion 

We have considered a class of pair-interaction stochastic processes in an A^-particle system 
and their associated pair interaction driven master equations. We have proved a chaos 
propagation theorem for this class of master equations and have used this result to study 
the kinetic limits of two biological swarm models, the BDG and CL processes. By inves- 
tigating the invariant density of the CL process, we have shown that the chaos property 
may be lost at large times. This work shows that the chaos property may be true even for 
processes that seemingly build-up correlations but may not be uniformly valid in time. 
Correlation build-up manifests itself at large time scales. In order to restore the validity 
of kinetic theory at these large scales, new theories must must be developed. This is a 
fascinating and widely open area of research. 
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